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We present a comprehensive study of primordial fluctuations generated from G-inflation, in which 
the inflaton Lagrangian is of the form K((f>,X) — G{4>, X)D<j> with X = — (<90) 2 /2. The Lagrangian 
still gives rise to second-order gravitational and scalar field equations, and thus offers a more generic 
class of single-field inflation than ever studied, with a richer phenomenology. We compute the power 
' spectrum and the bispectrum, and clarify how the non-Gaussian amplitude depends upon parameters 

s ! , such as the sound speed. In so doing we try to keep as great generality as possible, allowing for non 

slow-roll and deviation from the exact scale-invariance. 
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I. INTRODUCTION 



Cosmological inflation [l| is now a widely accepted paradigm explaining the flatness, homogeneity, and isotropy of 
the observed Universe. In the most common scenario, inflation occurs when the inflaton, a scalar field driving the 
accelerated expansion, rolls down a nearly flat potential slowly. During this slow-roll stage fluctuations in the inflaton 
field are generated quantum-mcchanically and stretched outside the Hubble horizon, which eventually reenter the 
Hubble radius in a later epoch to be a seed for the large-scale structure of the Universe. The detailed shape of the 
potential can be probed by observing the power spectrum of fluctuations in terms of the cosmic microwave background 
(CMB) anisotropies 0- As to theoretical approaches, much effort has been made to determine the inflaton potential 
in the particle physics context. However, single-field inflation with a canonical kinetic term and a nearly flat potential 
' is not the only option to induce the accelerated expansion and to produce almost scale-invariant perturbations with 
an appropriate amplitude. Liberating inflation models from the standard assumption, one may consider a variety 
of interesting scenarios: multiple scalar fields might participate the inflationary dynamics, the kinetic term of the 



inflaton(s) might be non-canonical p^, and a scalar field other than the inflaton might be responsible for the density 
perturbation From a high-energy physics point of view, supersymmetric theories naturally provide many scalar 
fields with flat potentials Q, and the Dirac-Born-Infeld (DBI)-type non-canonical kinetic term naturally arises from 
^ . D3-brane motion in a warped compactification Q. 
. i-H ' Different inflationary scenarios can be distinguished by future and on-going experiments such as Planck [7j , aiming 
, to obtain better constraints on the amount of non-Gaussianities in the primordial curvature perturbations as well 
H ' as on the spectral index n s , its running, and the tensor-to-scalar ratio r. The standard canonical slow-roll inflation 
models produce negligible non-Gaussianity 0] , while exotic inflationary scenarios are expected to predict measurable 
non-Gaussian signals. In the context of single-field inflation, non-Gaussian perturbations have been computed for the 
Lagrangian of the form 0, [HI 

Ci=K(<f>,X), (1) 

where <fi is the inflaton and X := — d ^4>d^ <f> / 2 . This class of models yields a sound speed c s different from the speed 
of light in general, and large non-Gaussianity is generated for c s -C 1. A significant non- Gaussian signal together 
with the confirmation of the consistency relation r = — 8c s 71t, where nx is the spectral index of primordial tensor 
perturbations, is a smoking gun of the inflaton Lagrangian ([I}. 
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In this paper, we consider a more general Lagrangian [TT], Il2| 

C 4> =K(<f>,X)-G{<t>,X)U<t>, (2) 

where K and G are some generic functions of the inflaton <j> and X. The new term G((f>, X)n<j> in the Lagrangian (0 is 
inspired by the Galileon interaction [13, [l^ an d reduces to the one having the Galilean shift symmetry, d^(j) ~^ d^ifi+b^, 
in the Minkowski background in the case G oc X. One of the most important properties of the Galileon Lagrangian 
is that the field equations do not contain derivatives higher than two. The interaction G(cj>, X)D<p is a generalization 
of the Galileon term XU<fs while maintaining the second-order property. In this sense, the Lagrangian @ defines 
a more generic class of single-field inflation than ever studied. Here, the Galilean shift symmetry is abandoned in 
exchange for generality, but one should note that the symmetry does not make sense already upon covariantization 
for any interaction that is Galilean invariant in the flat background. 1 (The name "Galileon" is therefore no longer 
appropriate when covariantized.) Cosmological applications of the Galileon interaction can be found in flBj with 
emphasis on dark energy and modified gravity. Primordial inflation based on the generic Lagrangian @ was first 
proposed very recently by [l2l . Il6j , and is dubbed G-inflation. Almost simultaneously the same Lagrangian was used 
to explain the late-time cosmic acceleration rather than the primordial one [H Ej. In the effective-field- 

theory approach (20l | was employed to see the consequences of imposing the approximate Galilean shift symmetry on 
the Lagrangian of primordial perturbations. Interestingly, the scalar field theory with the term can violate the 
null energy condition stably. This fact motivates the authors of Refs. [2l|, [22[ to propose a radical scenario of the 
earliest Universe alternative to inflation. Some specific form of the above type of interaction arises from a probe brane 
action in higher dimensions [23| and from the Kaluza-Klein reduction of Lovelock gravity [U, [25| . A supersymmetric 



completion of Galileons is explored in [26 1. 



The purpose of the present paper is to understand the nature of cosmological perturbations generated from G- 
inflation. We rederive the power spectrum and the tilt of the spectrum without assuming slow-roll, clarifying how 
the (approximate) scale-invariance is achieved in G-inflation. We then calculate the cubic action for the curvature 
perturbation and evaluate the full non-Gaussian amplitude, again without assuming slow-roll and the exact scale- 
invariance. Throughout the paper we try to make our formulas as general as possible, which we hope maximizes the 
usefulness of the results. Recently, non-Gaussianity from G-inflation was calculated neglecting a number of terms 
working in the de Sitter limit [27j and in the slow- roll limit [28j . See also a recent work by Naruko and Sasaki, in 
which the superhorizon evolution of the nonlinear curvature perturbation from G-inflation is addressed (29| . 

This paper is organized as follows. In the next section we review the basic properties of G-inflation and derive the 
power spectrum of the curvature perturbation. In Sec. Ill we compute the cubic action for the curvature perturbation 
to evaluate the three-point function in G-inflation. 



II. G-INFLATION 



We start with a brief review on the basics of G-inflation [12, [lj| . The scalar field Lagrangian for G-inflation is given 
by Eq. @. Assuming that cj> is minimally coupled to gravity, the total action we are going to study is 



S 



d 4 Xy 



-fJ 



M 2 



In the following we will set Mp\ = 1. The energy- momentum tensor T^ v of the scalar field is given by 
T„ v = K x Vp<j& v <j> + Kg^ - 2V (|i GV„)^ + g^V \GV X (p - G X ^^A- 



(3) 



(4) 



Here and hereafter we use the notation Kx for dK/dX etc. Varying the action with respect to (f>, we obtain the 
scalar field equation of motion, 



K X U4> - X X x(V M V^)(V^0VV) - 2K 4>X X + - 2 (G - G+xX) U<j> 
+G X [(V^V^XV W) - (n^) 2 + fl M „WV>] + 2G 0A -(V M V^)(V^V^ 
-G xx (VV^ - g^Ucfi (V M V0) V^V A ^ = 0, 



'i.GssX 



(5) 



Concerning this point, one may worry about the naturalness of G-inflation models discussed in the present paper because there is no 
symmetry to protect the Lagrangian. However, it should be noted that symmetry, if present, must be broken at least to end inflation. We 
therefore will not provide a symmetry-based argument but rather take a phenomenological approach, assuming that some UV complete 
theory would give the (in some sense fine-tuned) Lagrangian that leads to second-order field equations. 
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which is of course equivalent to the conservation equation Vj/T " = 0. One verifies from Eqs. ((4]) and ([5]) that the 
gravitational and scalar field equations are indeed of second order. 

Higher order Galilcon terms (with a (^-dependent coefficient) such as f(<p)X [2(D0) 2 — 2V^V ! ,0V M V ly (/> + RX\ can 
be added to the scalar field Lagrangian while keeping the field equations of second order. Although the effect of such 
higher order Galileons might be interesting in the context of primordial inflation, we leave the issue for future study 
and concentrate on the Lagrangian of the form ([2]) in the present paper. 

A. The background equations 

Let us consider homogeneous and isotropic background: 

d.s 2 = -dt 2 + a 2 (t)dx 2 , (f> = <j>(t). (6) 

Although the energy-momentum tensor (jj) cannot be recast in a perfect-fluid form in general (Tlj . for the above 
cosmological ansatz it takes the desirable form = diag(— p,p,p,p) with 

p = 2K X X - K + 3HGx<p 3 -2G<t,X, (7) 

P = K-2{g<i, + Gx'<£)X. (8) 

The gravitational field equations are thus 

3H 2 = p, (9) 

-3iJ 2 -2H = p, (10) 
and the scalar field equation of motion is given by 

K x (j> + ZH^j + 2K XX X^ + 2Kx4,X - K$ - 2 (G* - G X4> X) + ZH^j 

+6G X [(HX)-+3H 2 X] - 4G X4 ,X4> - 2G 4>4> X + 6HG XX XX = 0. (11) 

If, for example, K is given by the standard, canonical kinetic term with a potential, K = X — V((f>), one can consider 
an inflationary scenario in which the energy density is dominated by the potential as in the standard case, while the 
dynamics of the scalar field is modified by the Gn<j) term, changing the potential that <f> effectively feels. This is the 
scenario proposed in [l6| and called potential driven G-inflation. Another possible scenario is that inflation is driven 
by </>'s kinetic energy which is kept almost constant with nontrivial functional form of K and G. In models with the 
exact shift symmetry, <f> — > (f> + c, i.e., K = K(X) and G = G(X), it is easy to obtain an exactly de Sitter background 
satisfying H = const and <f> = const. This may be regarded as a generalization of k-inflation Q, and we call the class 
of models kinematically driven G-inflation [12J . Deferring the summary of these two specific classes of G-inflation to 
Sec. Ill Cl we now move on to describe the general properties of the power spectrum of primordial perturbations from 
G-inflation. 

B. Power spectrum 

In this section we derive a series of general formulas for linear cosmological perturbations without assuming any 
specific form of K and G. Wc work in the unitary gauge, 4>{t, x) = 4>(t). 2 Using the remaining gauge degree of freedom 
the linearly perturbed metric is taken to be 

d.s 2 = -(1 + 2«i)di 2 + 2a 2 d i Pidtdx i + a 2 (l + 2"K)dx 2 . (12) 

Expanding the action to second order in perturbations and then varying with respect to a\ and we obtain the 
following constraint equations: 

Tl = 6ai, (13) 

d 2 

— in + a 2 Qp 1 ) = XQa u (14) 



The unitary gauge does not coincide with the comoving gauge, STf = 0, in the case of G-inflation [T3 |. This fact stems from the 
imperfect-fluid nature of the energy-momentum tensor ((4J. 
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where d 2 := S tj didj, 

6 := H-$XGx, (15) 
g := K x + 2XKxx + 6G x H<j) + 6G x X 2 -2(G 4> + XG^ x )+6GxxHXj). (16) 

Substituting the constraints (JT3J) and (|T4")) to the action, we arrive at the quadratic action for 7Z [ll], [l2| : 



5 2 = / dtd 3 x a 3 a 

where 



_1 

a 



(17) 



c 2 s ■= f , (18) 



and 



(19) 



T := K x + 2G A - [4> + 2H^ - 2G X X Z + 2G X xX4> - 2 (G* - XG^) . (20) 

One can verify that setting G(4>, X) = the quadratic action (jTTJ) reproduces the expression obtained for It- 
inflation [3(| ■ It is useful to notice that a can also be expressed as 



e 4>xg x 
e 2 e 

Let us define three parameters that characterize the rate of change of three background quantities 



a = -^2+ o ■ ( 21 ) 



In this paper we assume that 



IP' S -=H^> 5 -=H^- (22) 



-k^ Ws^ m*°> ™ 

but we do not neglect e, s, and 5. (In the next section, however, we will assume some stronger conditions to evaluate 
the bispectrum.) It should be noted in particular that a is not necessarily small, in contrast to the usual (k-)inflation 
models in which a is degenerate, i.e., a = e < 1 [33]. Even in the slow- roll limit we may have a > 1 in G-inflation. 

Under the assumption that the parameters defined in (|22[) are constant (but not necessarily very small), it is 
straightforward to solve the equation of motion derived from the action (|17l) and compute the power spectrum of 
1Z fl2j . For this purpose it is convenient to define a new time coordinate y by dy = c s dt/a (3lj . In terms of y, the 
scale factor, the sound speed, and a are written as 

c s*{y/y*) 1 ^ 1 c s ' , , N_ s /(i_ e _ s ) , , ,-s/n-e-s) /OA \ 

a = i — \> c B = c s *y/y„) ; , <T = <j«{y y*) n > ( 24 ) 

{-y*)H*(l - e - s) 

where t he qua ntities with * are those evaluated at some reference time y ~ y*. Using a new variable u := zlZ with 
z := a\j2o J 'c s , the equation of motion can be written in the Fourier space as 

~z" 



z 



< + [V - — )u k = 0, (25) 



where the prime denotes differentiation with respect to y and we find 

- ( a/2-, 2" g 2 - 1/4 3 - e - 2s + j 

z oc -y 7 9 , — = - , with q := — — . (26) 

z y A 2(1 — e — s) 

The normalized mode solution to Eq. (|25j) corresponding to the Minkowski vacuum in the high frequency limit is then 
given in terms of the Hankel function by 

R ■ r(i) 



uk = ^V=yHW(-ky). (27) 



5 



Wc thus write the operator 1Z using the creation and annihilation modes as 

fc(k,y) = ip(k,y)a- k + ip*(~k,y)a f _ k , 

my) - Uk[y) 



with the commutation relation [dk, at,] = (27r) 3 (^ 3 )(k — k'). This immediately leads to the power spectrum [l2j], 



R 



2^ 



2 2q 



-3 



r(«) 



r(3/2) 



(1 - e - s) 2 H 2 



The scalar spectral index is found to be 



1 = 3 - 2q 



4tt 2 



2e + s + 5 
1-e-s ' 



2(TC S 



fey=-l 



(28) 
(29) 

(30) 
(31) 



The above formula has been derived without assuming the smallness of e, s and S, though we have assumed that they 
are constant. In this sense, the above expression is more general than that given in [lj, [l(| [2?], [H|. To ensure the 
scale invariance we require 2e + s + 5 ~ 0. However, this does not force each parameter to be as small as 0{n s — 1); 
each can be large, e,s,(5> 0(n s — 1), but the three may cancel each other out to produce an almost scale- invariant 
spectrum. This possibility was first pointed out by [3l[ in the less generic context of DBI inflation, for which a = e and 
consequently 6 = 0. We leave this interesting possibility open, and will complete the following calculation without 
taking the slow-roll limit. We would stress again that even if we consider the slow-roll limit, a is not necessarily 
slow-roll suppressed. 

Since the inflaton field is minimally coupled to gravity, the nature of tensor perturbations is the same as the standard 
one and is dependent only on the geometrical quantity H = H(t). In the slow-roll limit, e = 0, the tensor power 
spectrum is given by Vh = 8(H/2ir) 2 . The tensor-to-scalar ratio r is thus given by 



r = 16cr c s , 

where just for simplicity the scalar power spectrum is evaluated also in the slow-roll limit, e = s = S = 0. 
For later convenience we introduce the following quantity: 



(j>XG 



x 



H 



or, equivalently, O = H(l — v). From Eq. (f2~Tj) we obtain 

v 



H(l-is) 2 1 



+ 



1 



For e = const, s = const, and 6 = const, the above equation can be integrated to yield 

1 a (y) , ( 1 1 a * \ („./.. ^(l+e)/(l- 



H 

e 



(y/y* 



l-v(y) 1 + e l + e + 6 ' 1 + e l + e + S / 

If we assume v = const then we have a = const. In this case the two quantities are related as 



1 



(32) 



(33) 



(34) 



(35) 



(36) 



Note in passing that the opposite is not in general true: for a = const Eq. (|34[) still admits time-dependent 



C. G-inflation examples 

1. Kinematically driven G-inflation 

Inflation can be driven by kinetic energy of 4>. This possibility was explored in (l2l |. Let us consider for simplicity 
the Lagrangian with exact shift symmetry <j) — > </> + c, i.e., 



K = K(X), G = G(X), 



(37) 
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and look for an exact de Sitter background satisfying H = const and = const. It follows from the field equations 
that 

3H 2 = -K, (38) 

K x + 3G x H(j) = 0. (39) 

For this background we have 

T = -^"(l-"), ( 4 °) 

„ K I XKxx XGxx\ , s 

g = -H i+ *- 2 ^f +2 ^fv< (4i) 

(42) 



l-i/' 



where ^ = <j>XGx/H = XKx/K = const. In evaluating the above equations we used the background equations (f38|) 
and (p]). 

The concrete toy model presented in [l2| is given by 

where M and y, are parameters. In this case, c s and a can be expressed in terms of /z. It turns out that the tensor- 
to-scalar ratio r = 16ctc s = 16ac s (a) is an increasing function of a, and a ~ v <C 1 is required in order for r not to 
exceed the observationally allowed value. Explicitly, one finds r ~ (8/v / 3)cr 3 / 2 ~ (16v / 6/3)(v / 3m) 3 ^ 2 Ell- 
Note, however, that v <§; 1 is not necessary to get a stable, prolonged de Sitter phase. As already emphasized 
above, a > 1 is made possible by a suitable choice of K (X) and G(X), provided that r = 16ac s remains not too large. 
In [27[ Mizuno and Koyama have studied the case with a ~ v <C 1 focusing their attention on the model (|4"3"|) . In 
contrast, the analysis in the present paper can apply to more general cases with a > 1. 

In the presence of exact shift symmetry the exact de Sitter solution is an attractor. Along this attractor the scalar 
fluctuations acquire an exactly scale- invariant spectrum. Making K and/or G weakly dependent on <j), one obtains a 
quasi-de Sitter attractor and thereby the spectrum can be tilted. Though we do not provide corresponding concrete 
examples here, more generic, possibly complicated, choices of K(4>, X) and G(</>, X) would lead to the interesting 
situation mentioned above: n s - 1 « 1 with e,s,i5 > 0(n s — 1). 



2. Potential driven G-inflation 

In fl6| a novel class of inflation models was proposed in which the energy density is dominated by 0's potential but 
its dynamics is nontrivial due to the GOcjj term. In particular, it was shown that slow-roll inflation can proceed even 
if the potential is too steep to support standard slow-roll inflation. The model examined in [l6[ is described by 

K = X-V{4>), G = -g{<f>)X. (44) 

For gV,), 1, the effect of the Gn</> term dominates in the slow-roll equation of motion for <j> and the potential is 
effectively flattened, leading to slow-roll G-inflation. In this regime one finds 

4 2 

a ~ -e and c s 2 ~ -. (45) 

Though a could be free from the slow-roll constraint in principle, in the present case it is actually related to e in a 
way different from standard slow-roll inflation. Since c 2 ~ const, the scale-invariant spectrum requires that e< 1, 
and hence a <C 1. 

III. BISPECTRUM 

In order to evaluate the bispectrum, we compute the cubic action for 1Z working in the ADM formalism (8l-[l0j. 

d.s 2 = - N' 2 dt 2 + h i:j (dx l + N l dt) (dx j + N j dt) , (46) 
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where 



a z (t)e 2n 5 i:j , N = l + ai+a 2 + 



Ni = a 2 d l {px+ & + ■■■) + N u + 



(47) 



with d l N n i = 0. Here, a„ and (3 n are (D(lZ n ). The fluctuation of the scalar field vanishes in this gauge. At linear 
order the above metric reduces to Eq. (JT3J). 

As pointed out in 0] , we only need to consider first-order perturbations in TV and N l to get the cubic action. (This 
holds true even in the presence of the GO(f> term.) Therefore, it suffices to use the first-order solution of the constraint 
equations, Eqs. (fT3)) and (| 14[) . supplemented with a vanishing first-order vector perturbation, Nu = 0. 

We plug the solution for a\ and /3i into the action and expand it to third order in 1Z. After cumbersome multiple 
integrations by parts, one ends up with 



5, = 



dtd 3 x a 3 



^ii 3 + c 2 nn 2 ^ 

H a i H z 



Cs d 2 n(dn) 2 



4 T K 2 d 2 K + C 5 HK 2 K 



i 2 H 2 



%d 2 n(dn ■ dx) + - 4 d 2 n(d x ) 2 + %n{dn) 2 

an a n 



(48) 



where \ := d 2 A with 



a := ^xgn = ^-n. 



(49) 



The dimensionless coefficients are given by 



Ci 

Co 

c, 
c A 
c 5 

c 6 
c 7 

C s 

Co 



H a 



a 
73 



1 + 2^ 



H 2 

c 2 e 2 



2<j>X (G x + XGxx) 



Ha H 2 a 



26 



H 2 cf)XGx 
63 ' 
2H 2 ^X{G x +XGxx) 
6 3 
H 2 5 

lie 2 



2c 2 H dt 
2H<f>XG x 

a 4>XGx 
4 " 



— 1 - e - 2s - 



a 
73" 



e 

IP a 

e 2 c 2 v 

2H a 
1T+2 



c 2 e 2 d& 2: 



26 \ 

HQ ' 



(50) 
(51) 

(52) 
(53) 
(54) 

(55) 
(56) 
(57) 

(58) 



where 



2 := 



XK XX 
2X 



2X< 



Kxxx + H4>G X + 6X 2 G 2 X + 5H<j)XGxx + 6X 3 G x G X x + 2H^X 2 G X xx 



(2G c j ) x + XG^xx) ■ 



(59) 



The last term is the field equation which follows from the quadratic action, 



5L_ 
SK 



dt 



(60) 
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multiplied by 



TV 



H 

cle 2 



4a 2 6 ; 



[~(dK) 2 + d^d^^diTidjTZ) 



2a 2 9 



- [d x ■ dTZ - d- 2 d l d j (diTidjx)] -(61) 



In deriving the above cubic action we have not performed any slow-roll expansion, so that we have kept full generality 
up to here. Taking the limit G — > 0, — > H, and a — > e, we can verify that the above equations reproduce the 
previous result derived for generic k-inflation models, = K((f> 7 X) [9l Hoj|. In particular, the C3, C4, and Cq terms 
are absent in that case. The C5 term is clearly a higher order term so that we will neglect it in the following. 
Employing the in-in formalism, the 3-point function can be computed from the following formula: 



d*'([ft(ki, t)K(k 2 , t)K(k 3 , t), H int (t')}), 



(62) 



to 



where to is some early time when the fluctuation is well inside the horizon, t is a time several e-foldings after the 
horizon exit, and the interaction Hamiltonian is given by 



#i„t (*) 



d xa 



H 



TV + C 2 Wl 



J2 



We use Eqs. (J2ZI) and 



to evaluate each contribution, which can be conventionally expressed as 

A 



(TZ^TZ^TZ^) 

A 



5> 



t 2+ k 3)^T3T3T3 ; 



HI- 



(63) 

(64) 
(65) 



M 



The power spectrum V-r. here is to be calculated for the mode with kt = fci + k 2 + ^'3- 

To proceed, we assume that v = const, which holds in a wide class of G-inflation models as described in Sec. Ill Cl 
We then immediately see that a = const, and C3, C§, and C7 are all constant in time as well. The coefficients are 
explicitly given by 



(I +<??(<? -e) 
(1 + 6)3 ■ 



C fi = 



2(l + g)((j-e) 



C 7 



4e - cr(3 - e) 
4(1 Tej 



(66) 



In order to evaluate the contributions from the TZ 3 (C\) and TZ 2 d 2 TZ (C4) terms, we further assume that I/Q and 
4>X 2 Gxx I H are of the form 



(j)X 2 G 



Q 

xx 



H 



Ql + -o) 



(67) 
(68) 



where J\,Ji-,Q\i and g 2 are constants. In kinematically driven G-inflation [12[ we indeed have I/Q = const and 
4>X 2 Gxx I H = const in the de Sitter limit. In potential driven G-inflation [l6| I/Q ~ const and Gxx = 0. Therefore, 
the assumptions made here are sufficiently general and reasonable. It then follows that C\ and C4 take the form 



c Vl 4 



C4 = ^4 + 4 



I ' 

2 



where 2?i , E\ , 2?4 , and £4 are constant and are given by 



£1 
84 



1 + e 
ajl + a) 

1 + e 
(1 + af 



1 + 2 Ji + 2 



(7 — £ + (1 4- (t)^! 



2J2 



(1 + e! 
9 (lW 



l + cr 



1 + CT 



£>1 



(1 - 2g 2 ) 



(69) 
(70) 

(71) 
(72) 
(73) 
(74) 
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Each contribution can now be evaluated as 
3 



A l 

A 2 
A, 
A 4 
A 6 
A 7 
A 8 
A 9 



2a 

1 

4 



(1-e-s) 



r(g) 



r(g) 


2 / fcift 2 ft 3 \" s 1 


r(3/2) 





r(3/2) 
1 c 3 



k 1 k 2 k 3 



2kf 



V 1 I l (n,-l) + ^-h(g') 



it ( n 3-e/l + a 
3I 2 (n s - 1) - -J- ( y— : 



(l-e-s)* 



r(«) 



r(3/2) 



2fc| 



W), 



-(l-e- S ) 

(7 



r(?) 



r(3/2) 
r(g) 



fcifc 2 fcs 
2*? 



P4/4K - 1) + %/4(g') 



C7 <T 



r(g) 



r(3/2) 
r(g) 



r(3/2) 



kik 2 k 3 



/ fcifc 2 fc 3 X 
2fc3 



4(9') 



2ft? 



W), 



r(3/2) 



2ft? 



Cg* 


r(«) 


2 


f fcifc 2 ft 3 \ 


~8~ 


r(3/2) 




K 2k? J 



-I 8 (n, - 1) + 



W), 



1 + e - 2s / 1 + a 



1 + e 



where c s * and Cg* are evaluated at sound horizon crossing, fc t y = —1, and 

s - 2e 



1-e-s 



The fc-dependent functions Im are given by 
k\klk\ (irz\ T(3 + z) 



h{z 

h(z 
h(z 
h(z 

Hz 
Hz 

Hz 



ft? 



cos 



(?) 



(?) 

2 + z 



i>j 
"KZ\ 2 + z 



(ki ■ k 2 )fc^ 

fc 2 ftfftf /ttz^ (6 + z)r(3 + z) 



(?) 



1 i*i 

r(i + z) + r(2 + z) 



T(l + z), 

kik 2 + k 2 k 3 + k 3 ki 



k\ 



(3 + z) 



k\k 2 k z 
k? 



ft? 



(ki • k 2 )ft : 



cos 



cos 



ft* 

(ki ■ k 2 )fc 3 2 
ft* 

-(f) (e 

(kl • k 2 )fc; 



(¥) 
(¥) 
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(3 + z)r(l + z) + (3 + z)r(2 + z)-i - r(3 + z)-% 

k t kf 



sym., 



r(l + z) + T(2 + z)- 



sym., 



ft/ 



1 — z k 



1 + 2. 



* 2 



-kik 2 k 3 



r(i + z), 



cos 



(3 + z)r(i + z) -r(2 + z) 



ft* 



sym.. 



In computing the above we have used the approximation 

i>(k, y) s y27r^ /2 (fc t )fcr 3/2 ^' 9 (l + e" lfca , |fty| « l. 



(75) 
(76) 

(77) 
(78) 
(79) 
(80) 
(81) 
(82) 

(83) 

(84) 

(85) 

sym., (86) 
(87) 
(88) 
(89) 

(90) 

(91) 

(92) 



One can check that setting a — > e the above equations reproduce the result of [32| . The field redefinition 1Z — > TZ„ + 
f(!Z n ) gives rise to the non-Gaussian amplitude proportional to S, which can be ignored in the present approximation. 
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The shapes of ^3(1, fe, ksj/kzkz and -Ae(l, fei &3)/fc2&3 are very similar to the equilateral one, as was already 
pointed out by The contributions from 1Z 3 (C x ) and TZ 2 d 2 TZ (C 4 ) give the same momentum dependence because 
the two terms are essentially equivalent after using the first-order equation of motion. 

The size of the three-point correlation function is conventionally parameterized by /nl defined as 



/nl = 30 



■k2—k3 



kf 



which can be computed straightforwardly by evaluating the amplitude A at k\ 
for /nl is given by 



k? = fca 



(93) 

fct/3. The exact expression 



h 



NL 



30 



r(3/2) 



3(1 



s) 



2a 



3(1 



sf 



V 4 



6 



+ ^/ 2 eqm (n s - 1) + 



(1 



2acl 
1 - e - s 1 



1 

(1 + ^)V 

3 



(1 + ef 



&c 2 



1 



1 



-/r qui (n s - 1) + 



3a 2 



£46 + 9' 



1 



/TV) 



2(1 + e) 1 + e 
3 



(3 



(»« - 2) 



(94) 



where we have defined i^ qul (z) := k t 3 lM(z)\k 1 =k 2 =k 3 , i- e -, 



I^ ni (z) 


:= cos 


f"KZ 




:= cos 


I ' TtZ 

\Y 


I^(z) 


:= cos 


l"KZ 

\~2 




:= cos 


CKZ 

\~2 



1458 



81 



2916 



243 

equi/.^ — T c< i ni ( 



(95) 
(96) 
(97) 
(98) 



and used the fact that 7 4 ° qu » = (1 + z/6)I^ ul (z), I 7 cqm (z) = r 2 qul (z)/2, r s qu \z) = 3/ 6 oqm (z)/(l - z), and 7 9 cqm (z) = 
I2 (^z) to shorten the expression. The above generic formula is involved, but an order of estimate of /nl is found to 
be 



/i 



NL 



o 



O G 



G 



x 



O or 



a := max{l, cr}. 



(99) 



This is one of the main results of the present paper. 

The sound speed at horizon crossing can be written in terms of kt k t . Under our assumptions 

we see that /nl can be expressed as /nl = fi + .f-2/c 2 *, where /1 and f 2 depend on e, s,a,T>i t 4, and £1,4, but are 
independent of k t . Therefore, the wavenumber dependence of /nl appears only through c 2 t , so that the tilt «ng is 
given by 



«NG - 1 



2s 



fl + /2Cs* 1 - £ 



(100) 



If e, s <C 1 and the main contribution to /nl is due to a small sound speed, then we recover the result of [IoT |. 
jing — 1 — — 2s, even in the presence of the GD</> term. 

We close this section by illustrating several examples of non-Gaussian shapes .4(1, k 2 /ki, k^/ ki)(k 2 / fci) -1 (£3/ fci) -1 
in G-inflation. First, let us consider the de Sitter limit of shift-symmetric kinematically driven G-inflation. The 
functions K and G may be written as 



K 



-X + Cl X 2 + c 2 X 3 + 



G 



X 
TP 



d 2 X 2 



d*X 3 



(101) 



11 




FIG. 1: The non-Gaussian amplitude .4(1, fe/fei, k^/k 1 )(k2/ki)~ 1 (k3/ki)~ as a function of A^/fci and k^/k\ for kinematicallly 
driven G-infiation. The amplitude is normalized to unity at an equilateral configuration, fo/fci = ka/ki = 1. The parameters 
are given by a = 0.36, c s = 0.03, g = 1, and X/<5 = 1, so that r ~ 0.17. The size of non-Gaussianity is /nl — 210. 




FIG. 2: The non-Gaussian amplitude A(l, ki/kx, k$ I k\)(ki j k\)~ x {kz I k\)~ x as a function of fe/fci and k^jk\ for kinematicallly 
driven G-infiation. The amplitude is normalized to unity at an equilateral configuration, fe/fci = k^/ki = 1. The parameters 
are given by a = 0.1, c s = 0.1, g = 60, and X/<5 = 1. The size of non-Gaussianity is /nl — 204. 



where Ci, di, and M are arbitrary in principle. Given v = <fiXGx/H and p = <fiX 2 Gxx / H , the former is related to 
XKx/K through the background equations, which in turn fixes the value of a. The latter is related to c^, but since 
the expression for c 2 contains both of the second derivatives Kxx and Gxx, Cg can be chosen independently of g. 
Third derivatives Kxxx and Gxxx appear only in the function X. In summary, in the case of kinematically driven 
G-infiation, the non-Gaussian amplitude in the de Sitter limit is completely determined by the four parameters 

cr, Cs, g, (102) 
y 

The four parameters can be written in terms of Ci, di, and M, but in practice the expressions are quite involved. We 
plot in Figs. HHS] the shapes of non-Gaussianity for different parameters. 

Another example of the non-Gaussian shapes we explicitly compute is given by potential driven G-inflation (l6l | in 
the slow-roll approximation. In this case we have a ~ 4e/3, ~ 2/3, ~ H (1 — e/3), and 



x = - 9 H4> + Qg 2 x 2 + -X 9<p 



(103) 
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FIG. 3: The non-Gaussian amplitude .4(1, k^/ki, k^/k 1 )(k2/ki)~ (k^/ki)~ as a function of k^/ki and k^/k\ for kinematicallly 
driven G-infiation. The amplitude is normalized to unity at an equilateral configuration, fo/fci = kz/ki = 1. The parameters 
are given by a = 0.1, c s = 0.1, g = 1, and T/Q = 300. In this case the shape peaks in the folded configuration ki = 2/C2 = 2&3. 



The contributions relevant at leading order in slow-roll are: 

A 1 ^]l 1 (0), A 2 ^-h 2 (0), A 3 ~~I 3 (0), A 4 c±h 4 (0), A 8 ~^-I 8 (0). (104) 
4 8 16 2 16 

More explicitly, the non-Gaussian amplitude is given by 
A 



1 h-^ &2 ^3 3/2 ^ > 7 2 j 2 \ ^ / 2 ; 3 1 

- i —-l[k t l. k ' k '-T4^ k - k ') 



(105) 



which is plotted in Fig. |U Taking the equilateral limit, the size of non-Gaussianity is found to be /nl = 235/3888 
0.06. The above result is insensitive to the inflaton potential. 



IV. CONCLUSION 



In this paper, we have studied G-inflation, i.e., generic single-field inflation obtained from the Lagrangian @. We 
have revisited the power spectrum and the spectral index to clarify how the (approximate) scale-invariance can be 
achieved in this class of inflation models, and determined the possible non-Gaussian amplitude without assuming slow- 
roll and the exact scale-invariance. The nonlincarity parameter /nl in G-inflation can be summarized schematically 
as 

/nl = O (i^j + O (^^f^) + ° . ° ■= max {l, (106) 

It should be emphasized that we have in principle no dynamical constraints that require a to be very small. If the 
first term dominates and a > 1, then we have /nl ~ 0-2 /c 2 - Therefore, large /nl an d large r = 16ac s are compatible. 
The situation should be contrasted with the models without the GO<f> term, for which a = e follows, and hence large 
/nl(~ 1/ c s) implies small r. 
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FIG. 4: The non-Gaussian amplitude .4(1, fea/fci, k- i /ki){k2/ki)~ 1 {k3/ki)~ 1 as a function of k^/kx and k^/ki for potential 
driven G-infiation. The amplitude is normalized to unity at an equilateral configuration, fo/fci = kz/ki — 1. The size of 
non-Gaussianity is /nl — 235/3888 ~ 0.06. 
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